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We reanalyze the neutrino events from SN1987A in 1MB and Kamiokande-II (KII) detectors, and 
compare them with the expectations from simple theoretical models of the neutrino emission. In 
both detectors the angular distributions are peaked in the forward direction, and the average cosines 
are 2 sigma above the expected values. Furthermore, the average energy in KII is low if compared 
with the expectations; but, as we show, the assumption that a few (probably one) events at KII 
. . . have been caused by elastic scattering is not in contrast with the 'standard' picture of the collapse 

\l ' and yields a more satisfactory distributions in angle and (marginally) in energy. The observations 

, give useful information on the astrophysical parameters of the collapse: in our evaluations, the 

^ ■ mean energy of electron antineutrinos is {E) = 12 — 16 MeV, the total energy radiated around 

' (2 — 3) X 10^^ erg, and there is a hint for a relatively large radiation of non-electronic neutrino 

5-H species. These properties of the neutrino burst are not in disagreement with those suggested by 

the current theoretical paradigm, but the data leave wide space to non-standard pictures, especially 
■^T ' when neutrino oscillations are included. 
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■ I. MOTIVATIONS AND CONTEXT 

m ■ 

o : 

. The detection of neutrinos from SN1987A marked the beginning of (extra)galactic neutrino astronomy 
' (see also |3,l3,i3l for comprehensive reviews of SN1987A observations). The observations of Kamiokande have been 
Q^, mentioned and recognized in the 2002 Nobel prize for Physics. However, when one studies the data, one meets a 
I ' number of surprising, unexpected or even puzzling features. Let us recall which are the main ones. 
O (1) The angular distributions of the events seen at Kamiokandc-II (KII) and at Irvine-Michigan-Brookhaven (1MB) 
> ' are more forward-directed than expected, for instance, the average cosines of the polar angles are (cos^^^^) ^ 0.3 and 
(cos 01^1^) -0.5. 

• • . (2) Also, the energy distribution of these two detectors seems not to be perfectly in agreement. In particular, (i?™) is 
ti^ilf of {E^^) (about 30 MeV), that is, a very marked difference even taking into account the different performances 
of the detectors. 

(3) Even the time distribution of the events in the two detectors looks to be different. However, when data are 
combined the distribution in time does not contradict the current picture of a 'delayed explosion' according to Lamb 
and Loredo analysis [To|. 

(4) The Mont Blanc events P occurred 4.5 hours before the other ones. This led some Authors to consider two-stage 
scenarios for the collapse 

In this work, we will focus on the discussion of the first issue and will stress the connections with the second one. 
More in general, we believe that these data raise several important questions that deserve attention, for instance: 
How likely is that the anomalies in the distributions are due to fluctuations, and in particular, how significant is the 
hint for some feature in the angular distributions? What can we learn (and what we can exclude) on the nature and 
the properties of the stellar collapse from these observations? 

A number of recent facts, beside the general considerations exposed above, testify the interest in having a fresh look 
at the SN1987A data: (a) several experimental evidences (in particular [i3,llJ,ll3) strongly suggest that SN1987A 
neutrinos oscillated in flavor; (b) the expectations of the emitted neutrino radiation has been recently reconsidered 



1Marialaura.Costantini@aquila.inf'n.itl 

fAldo.Ianni@lngs.inf'n.it" 

• |Francesco.Vissani@lngs.infn.it| 



2 



[16)|, suggesting a new paradigm for the distribution of neutrinos and antineutrinos; (c) there have been improvements 
in the description of the cross section of 

De p ^ n, IBD reaction (from 'inverse beta decay') (1) 

in the energy range relevant for supernova neutrinos |l7tll8l |. Moreover, it is correct to recall that we do not understand 
yet the theory of core collapse supernovae, and therefore one could argue that we miss the most important ingredient 
for a proper interpretation. However, a reasonable working hypothesis is to describe the emitted neutrino radiation 
by a model with few parameters, suggested by the 'delayed explosion' scenario proposed in 0, see |23| for a recent 
report. This is the point of view we will adopt in a large part of the present investigation. 

We will describe and motivate in the rest of this Section what we assume (based on expectations and observations) 
as a reference neutrino flux. We will discuss a standard (but updated) comparison of observations and expectations 
in Section 2, based on IBD hypothesis (see below), that will permit to further define the parameters of the model 
of neutrino emission. In Section 3, we will use this model to analyze the angular features of the spectra, state the 
situation quantitatively, and consider a few alternatives to improve the agreement with the data. We summarize the 
results obtained in the last Section. 



A. Neutrino flux 

A simple model of the fluxes [T^ of supernova neutrinos attributes the following spectra (with three different average 
energies Eq) to any species ly^, or v,j.-x being any among muon and tau (anti)neutrinos: 

F N 

= ^ ;^ ^"e-("+l)^ z = ElE^ (2) 

where z = e, e, a; and iV = (a + 1)"+Vr(a + 1). The total fluence at the detector is / E^i(E)dE = Ei/AiTD'^, thus £i 
is the amount of irradiated energy in the neutrino species i (the flux is supposed to be emitted isotropically) and D the 
SN-detector distance. Numerical calculations find that the time integrated flux usually called 'fluence', is rather 
well described by this ansatz; in particular, the deviations from a thermal shape are not large, and can be described 
as we do here by setting a = 3 for all neutrino species (a = 2 amounts to a Maxwell-Boltzmann distribution) . Finally, 
the meaning of Eq is just the average energy of the species considered (we will take Eq — (Eg) in the following). 

The total energy emitted in neutrinos can be estimated by simple considerations. In fact, the total gravitational 
energy irradiated is £b ^ •^GNM^^/bRns, and using for the neutron star a mass of A/„s = (1 — 2)Mq and a radius of 
Rns — 20 km (Mo/M„s)^/'^, we get £b ~ 5) x 10^'^ erg. The amount of energy that goes in the specific flavors 
is uncertain. Since Sf. is not very important for the observed signal (see below), we will always set £e = fg, unless 
stated otherwise. Instead, we will distinguish three cases for the emitted energy £x (assumed to be equal for u^, i/^, 
Vr and i>f, so that £b = £e + £e + '^£x) that, as we will see, plays a more important role: 

1. £x — £e- This is the so-called 'equipartition', often adopted in theoretical analyses. 

2. £x = i?e/2: This is the case when a large part of the radiation goes in electron neutrinos. 

3. £x = 2£g: Finally, in this case most of the radiation goes in muon or tau neutrinos. 

The average energies are important parameters. {E^) is greater than (Eg), but the amount of hierarchy found 
in modern calculations is not very large. A typical ratio is in the range 1 — 1.2. In the following, we will assume 
(unless stated otherwise) (E^) — l.l(£'g). The average energy of the electron neutrinos (Eg) instead is not of crucial 
importance for the observed signal. It can be evaluated by prescribing a condition on the emitted lepton number 
ALe = N{Vf,) — N{i?e) where N{ve) = £&/ {Ee) and similarly for the antineutrinos; we will assume that the electrons 
contained in one solar mass of iron are converted in neutrinos. 

The crucial parameters needed to describe the neutrino signal are the antineutrino average energy, 

Eq = {Ee) ==12-18 MeV (expected) (3) 

and the total energy irradiated in antineutrinos 

£e = {2- 10) X 10^2 erg (expected) (4) 

Both of them have considerably uncertainties, especially the second one. For this reason, the uncertainty in the 
distance of the supernova is usually considered unimportant; here, we will assume D — 52 kpc, and discuss this point 
later. 
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B. Impact of neutrino oscillations 



Motivated by the solar and atmospheric results, we assume that the three neutrinos i^e, i'^ and i^r have mass and 
mix among them. Following simple minded theoretical expectations, we will further assume in most of this paper that 
the heaviest state is separated by Ato|]^ « 2.5 x 10"'^ eV^ from the other two, whose splitting is Atoj^ ~ 7 x 10~^ eV^. 
The known mixing angles are 023 ~ 45°, 9i2 ~ 34°, while ^13 is unknown but presumably it is not very small (we 
take « 6° when needed, but its impact on the oscillations is usually of minor importance). With these parameters, 
the emitted fluxes from SN1987A, described in Sect II Al should be modified to account for the MSW effect in the 
star mmiHlillill (among first papers on the topic, we recall [H IH [s^] ) : 

\ $g ^ Pee*e + (1 - Pee)^. ^ ' 

where the two probabilities of survival are Pgg = cos^ 6*12 « 0.7 and Pee = sin^ 6*13 « 0. (We will discuss other 
possibilities, as a very small ^13, or an 'inverted' mass hierarchy, in Sect lmPl l The MSW effect of the Earth modifies 

Pee 

by a minor amount |3l| . 

Two remarks are in order: (1) It is difhcult to conceive that oscillations did not occur; for instance, the MSW effect 
related to solar Am^ happened unless there was a drastic modification of the mantle of the star for densities around 
10 gr/cc, which seems unlikely. (2) The effects for D^. are of about 30 %, while those for i/^ can be much larger; for 
instance, in the framework described above, the observed Ve flux corresponds almost exclusively to emitted or Vr- 
This is the reason why has an important role for the signal seen in terrestrial detectors. 



II. THE IBD HYPOTHESIS 



In the model previously described and with the expected values of the parameters, it is a fact that most of the 
events expected at KII and 1MB are due to the inverse beta decay process. This is the reason why several analyses 
adopted the simplifying hypothesis that all events come from IBD (see e.g. 112 )• We begin by repeating such a more- 
or-less standard analysis, with three specific aims: (1) stressing observables with a clear physical meaning (rather than 
attempting a global analysis of the data); (2) discussing how the various data fit in the theoretical picture; (3) getting 
more specific values of the parameters of neutrino emission. The calculations of the expectations are quite simple. 
In a detector with Np protons and with detection efficiency e{E) (function of the positron energy) one integrates the 
differential event rate 

over the allowed range, obtaining the value of the observable of interest-e.g., the visible energy in Cerenkov detectors 
Ems — E (note that the fluence $ should be thought as differential in the transversal surface and ui Ei,). In Figure ^ 
we show the value of two simple but interesting observables: the mean energy and the number of events. In this 
Figure, the effect of varying (E'g), {Ex) and 8x in certain ranges are illustrated. For 1MB, we reduced the expected 
number of events by 13% to take into account the dead-time occurred during the detection of the burst 0- Let us 
comment the results in some detail: 

a. Average visible energy This observable has the advantage of being independent of the total energy emitted, 
and of having relatively small errors: 



I {ED-{E,,s)^ _ r 2.4 MeV at KII , . 

"^"""V TV ~ \ 2.6 MeV at 1MB ^' > 

in the IBD hypothesis Eyis is the energy released by the positron. While 1MB points to a range of values nicely 
consistent with expectations, compare with eq.©, the data of KII point to somewhat lower values. Note that we 
discard from the analysis the sixth event of KII, since it has Ni^a — 16, below the threshold Ni^n = 20 of software 
analysis Q ■ Indeed, it should be remarked that the lower energy events are those for which pollution from background 
is more likely; in particular, in the window of 12 sec in which the supernova neutrinos have been detected at KII, we 
estimate an average of about 2 background events 0. From Figure [2 one sees that the impact of a variation of {Ex) 
and £x on the expectation for the average visible energy is not large. 
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TABLE L Approximate coejficients of the angular distributions for IBD using eg. J^) . 





fflo 


ai 


02 


KII 


0.499 


0.030 


0.003 


1MB 


0.495 


0.104 


0.014 


1MB with bias 


0.492 


0.154 


0.024 



b. Number of events The observables N =11 and TVi^t^ — 8 have large Poisson errors, but permit to estimate 
the energy emitted from the supernova (whereas, the previous observable is not useful for this purpose). In the two 
plots on the right of Figure the energy emitted in any species of neutrino is chosen by default to be 4 x 10^^ erg. 
We see that the agreement with the expectations is not bad, and also that the impact of a variation of is not fully 
neghgible. This is easy to understand and to keep into account; indeed, the signal scales roughly as 0.7£g + 0.3£x, 
thus a variation in S,^ can be well simulated by a variation of the total emitted energy. 

c. Summary Using these results as a guide, wc further specify the parameters of the model and assume 

Eo = {Es) = 14 MeV, = 4 x 10^^ erg (8) 

These values should be thought as compromises between contrasting needs. For instance, KII energy spectra suggest 
lower values of i^o, whereas in order to reproduce the number of events at central values, we need Eq ~ 18 MeV and 
f e = 3 X 10*^^ erg. In view of this situation, we find it reassuring that the values shown in eq.® are in accordance 
with the expectations from a 'standard' collapse (as defined in eqs. Q and Q). 

We conclude this first part of the analysis reassessing that within the 'standard' model of the collapse and with the 
parameters of eq. the observed average visible energy at KII looks a bit small (see again Figure ^1. 



III. WHAT IS THE MEANING OF THE FORWARD EVENTS? 



In this Section we study the angular distribution of the events from SN1987A. Thus, we select the events from the 
two water Cerenkov detectors operative at that time, and use the data from for 1MB, and those from Q for KII. 
Both angular distributions are rather forward-directed. To state this more precisely, we calculate the average angles: 
(cos^^^^) — 0.29 ± 0.27 and (cos^^'^^) — 0.48 ± 0.34. Here we have used a weighted average and the corresponding 
standard deviation errors. 

d. Beyond the IBD hypothesis We can compare the data with the expectations from the IBD hypothesis. Using 
parameterized angular distributions 

dN/dz — oq + aiz + 022;^, where z — cos9 and as in Tab^ (9) 

obtained from [l^ we find that both central values are above the expected ones: 2.3cr for 1MB and 1.7cr for KII. 
This conclusion is in agreement with ^17) . In this study, we adopt the model defined in the previous Section with the 
parameters of eq. We checked that a variation of these parameters is not crucial for the conclusions, while £x is 
of greater importance. It is simple to explain the reason: The only type of events that is strongly forward (and thus 
is able to affect the angular shape of the distribution) are those from 

Vi e ^ Vi e, ES reaction (from 'elastic scattering'), i = e,e,x (10) 

This reaction receives contributions from all neutrino types, and v^. gives the largest one. But due to oscillations, 
eq. ((Sj, the observed Ve flux is originally due to v^', this implies the relevance of £x, namely, the energy emitted in 
i^^.r- The hypothesis that one or more forward peaked elastic scattering events could be present in the data samples 
of 1MB and KII has been already considered in the past, see e.g. 7, 33j .Mi .2^ .2^ i22}- In this analysis, however, 
we update the angular distributions for IBD events and the model for neutrino emission, compare different statistical 
inferences and include oscillations with recently measured parameters. 

e. Instrumental effects A point to take into account is that the angular distributions (and in particular the one 
of ES) are modified in an important manner by instrumental effects. This is due to multiple scattering and limited 
angular resolution of the detectors, and it is called 'smearing' of the angular distributions. In order to account for 
this, we use the following distribution |33 |: 



Psmicos 9s)d cos 9s =iVe"^'/^'^' sinOgdes 



(11) 
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FIG. 1; Comparison of observations (horizontal strips) and expectations calculated in the IBD hypothesis. The left panels show 
the average visible energy, the right ones the number of events. The upper panels are for KII, the lower ones for 1MB. For 
any panel, we show 4 expectations curves. The continuos (red) ones correspond to a variation of the theoretical parameter Ex 
(—energy radiated in v^^r) in the range (2 — 6) x 10^^ erg. The dashed ones (green) correspond to a variation of the theoretical 
parameter (Ex) (=average energy of the emitted v^^, i/p, Vr and Vr) in the range (1 — 1.2) x _Eo. 



where ds is the angle of smearing and Og is a measure of the effect; N in eq. (|ll|l is just a normahzation factor. For 
KII, where the smearing is shghtly more important, we choose Ug in such a way that the mean angle B from Hll|) 
corresponds to the mean error 60 determined from the data '?|. For the whole set of data we have 59 ^ 25°, while 
considering only data with 9 < 30° we find 69 ^ 18°. As a consequence we decide to study the two cases, when 
(Ts = 15° and ~ 20°, for which 9 ~ 18° and ~ 25°, respectively. Based on similar considerations 0] we set as — 16° 
in 1MB. Using eq. Hll|l and p^s e.g. from 41] we can determine the reconstructed angular distribution as follows: 

Per(n • m) = / (fp psm{m ■ p) pes(n • p) (12) 
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FIG. 2: Reconstructed angular distribution for elastic scattering events m KII. 



where n, m and p are unitary vectors for the SN, the reconstructed and the emitted direction respectively, n-m = cos6', 
m • p = cos 0s, n • p = cos cos 0s + sin sin 0^ cos 0s , and dP'p — d cos 6 sd(j)s/4:TT is an element of solid angle from 
which the signal receives a contribution. The reconstructed distributions we obtain for KII are plotted in Figure |5J 
Since the angular distribution pes of ES is rather narrow (especially when taking into account detector efficiencies) 
the reconstructed angular distribution of ES is mostly dictated by instrumental effects. 

A. Angular distribution of 1MB 

The normalized positron angular distribution for inverse beta decay is usually taken in the simple approximation 
dN/d cos 9 — 0.5 + acos^; in particular, in the 1MB report it is assumed a ~ 0.07. In the same paper it is 
pointed out that to account for the experimental polar-angle efficiency, one can introduce a 10 % angular bias. This 
is equivalent to replace dN/dcos9 ^ (I + O.I cos9) dN/dcos9. We use the improved cross-section for IBD from 
to determine the parameters {i = 0, 1, 2) in Tab.^ that enter in the angular distribution of eq.@. We notice that 
ai's in Tab^do not depend significantly on the assumed mean energy Eq. 

We have used eq. © to test the hypothesis the data from 1MB come from IBD events, employing the Smirnov- 
Cramer-Von Mises (SCVM) statistics 39] . As shown in the left panel of Figure O the goodness of fit (g.o.f.) for this 
hypothesis is equal to 6.4%. The improved IBD angular distribution changes the previous result (4.5% 0) by only 
a small amount due to the poor statistics. However, the importance of using the improved angular distribution is 
evident when we compare the old significance without angular bias with the new one, since 1.5 % i increases to 
4.2 %. In the same Figure we show the cumulative distribution in the hypothesis of having 1 ES event in 1MB. 

To study the possibility to have a small contribution of Elastic Scattering (ES) events in 1MB and later in KII we 
have exploited the Maximum Likelihood (ML) method "4^. In this framework the likelihood function is written: 



L 



^obs 



= Jl P ( cos^^; 



^obs 



(13) 



where n/riobs is the parameter which measures the fraction of ES events, riobs being the total number of experi- 
mentally observed events for the SN. The angular distribution p{cos9i]n/nobs) can be written as p{cos9]n/nobs) — 
{n / riobs) PeT i'^'^^ ^) + (1 ~ n/nobs)pibd(<ios9), where p^f^ and pibd are the angular distributions for ES and IBD, re- 
spectively. It turns out that in 1MB the best-fit is found for n/uobs = 0. The effect of the smearing is not particularly 
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Scattering angle for SN 1987A in 1MB Scattering angle for SN 1987A in Kll 




-0.75 -0.5 -0.25 0.25 0.5 0.75 1 -1 -0.75 -0.5 -0.25 0.25 0.5 0.75 1 



COST? COST? 

FIG. 3: Theoretical and experimental angular distributions in 1MB (left panel) and KII (right panel). The 1 sigraa range and 
the goodness of fit figures are also shown. For KII, we adopt Os — 15°. 



important in 1MB. In order to determine an upper limit on the likelihood parameter we have built a posterior proba- 
bility distribution function (p.d.f.) by normalizing the likelihood function and considering a uniform prior p.d.f. which 
is equal to one for n/nobs > 0, zero elsewhere. It turns out that n/nobs < 0.12 at 68.3% C.L., namely the 1MB angular 
distribution admits one ES event at most. 



B. Angular distribution of KII 

As stated above we consider nobs=ll out of 12 candidate events 3 and assume the event number 6 due to back- 
ground. In Figure 121 we show the reconstructed angular distributions for the cases tJs = 15° and as = 20°. The 
smearing effect in KII plays an important role (without smearing, {0) = 10°). Using the ML from eq. H13|l we have 
computed the likelihood ratio L[n/nob.s) / Lmax with n = 0, 1, 2, 3 to quantify the probability to have zero, one or more 
ES events on the basis of the angular distribution. In Figure^we show the normalized ML function against n/riobs- 
Minimizing — \TiL{n/nobs) the best-fit is found for n/nobs = 0.35 ± 0.20(1(t) for as = 20° and n/nobs — 0.23j;g'^g(l(T) 
for as — 15° . So, the ML test suggests that a few ES events are present in the KII dataset. 

As for 1MB we have exploited the SCVM test (see Figure 13). Moreover, we have worked out the probability to 
have n = 0, 1, 2, 3 ES scattering events using the expectations based on the SN model described above. For this latter 
case we have written the probability to have n ES events out of a total of nobs—^^ as the product of two Poisson 
distributions, that is equal to: 

Pin) = 71,^;"= X Bp{n, nobs) (14) 

nobs ' 

where the first factor is a Poisson distribution with a mean value n^xp — n^s + nibd and the second one a Binomial 
distribution with a trial probability p = n^s/nexp ~ 0.03; for instance, for the equipartition scenario we found 11.9 
IBD events and 0.39 ES events. (Incidentally, one should notice that the calculation of the ES number of events is 
very sensitive to the experimental efficiency at the lowest measurable energies). In Tab.HTlwe summarize our results. 
Similar calculations were made for 1MB. 
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FIG. 4: Normalized likelihood function versus the fraction of expected ES events in KII. For the smearing of the angular 
distribution we show the cases a = 20° and a = 15° . 



C. Summary for the 'standard' scenario and remarks 

It is instructive to compare here the outcomes of the two statistical tests we used: the Smirnov-Cramer-Von Mises 
(SCVM) and the maximum Hkelihood (ML). The comparison with the data follows completely different strategies: 
the ML method is 'local' in the sense that it profits of events that fall under the ES bell (of FiglSJ, while the SCVM is 
'global' in the sense that it tries to minimize the maximal distance between the theoretical curve and the observed one. 
For 1MB the SVCM test suggests more elastic scattering events than the ML test does (the reason can be understood 
from the right-most part of figure the most forward event has polar angle 33° ± 15°, thus its central value is not 
forward enough to suggest an ES event). Instead, for KII the two tests give very similar indications. 

Next, we take into account also the theoretical expectation on the number of ES events, and use it together with 
the results from the angular distribution. We combined the information from the SN model and that from the SCVM 
analysis in Tab.HTI multiplying the probabilities and normalizing the resulting distribution to one. The results for KII 
are shown in Tab. IIIII and the combined probabilities are given for the case a = 20° (the case tr = 15° gives about 
the same result). In particular, from Tab. IIIII we see that one ES in KII can be accepted at about the same level we 



TABLE II: Expected probabilities that zero, one or more events in KII are due to ES, for as = 15° (upper part) and as = 
20° (lower part). The first line shows the a priori expectation from the model of eg. I^, while the second and the third line use 
the information from the observed angular distribution. 







n = 1 


n = 2 


n — 3 


from SN model 


50.5% 


38.9% 


9.7% 


0.9% 


g.o.f. from SCVM 


8.6% 


26.7% 


58.5% 


81.4% 


likelihood ratio L{n)/Lmax 


0.35 


0.73 


0.97 


0.98 


from SN model 


52.3% 


37.5% 


9.2% 


1.0% 


g.o.f. from SCVM 


8.6% 


24.9% 


53.8% 


87.6% 


likelihood ratio L(n)/Lm.ax 


0.14 


0.39 


0.69 


0.92 
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TABLE III: Relative percentage probabilities to have a given number of ES events in KII dataset, estimated from observed 
angular distribution and theoretical expectation on the fluxes, for 3 hypotheses on . 





n = 


n = 1 


n = 2 


n = 3 






52.3% 


37.5% 


9.2% 


1.0% 




= S-e 


40.0% 


42.3% 


15.3% 


2.5% 




= 2 


29.0% 


43.6% 


22.3% 


5.1% 



could accept zero events. Moreover, even the probability to have two ES events is indeed not negligible. Repeating 
the exercise for 1MB, we find that at 'equipartition' (fg — Ex) the combined probability to have zero (one) events is 
80 % (19.9 %). 

Some remarks are in order. 

(1) The data of KII show that the most directional events have energies above 20 MeV. Taking this experimental fact 
into account, we checked the probability to have events with Ef. > 20 MeV for the three scenarios of SN considered. 
As shown in Figure this probability is about 16%. So, it seems not unlikely |45| from this point of view to have 
measured ES events with energies above 20 MeV. 

(2) The presence of one or more ES events in KII dataset goes in the right direction to explain the disagreement 
between 1MB and KII average energies. However the effect is admittedly small, since for instance 1 ES event that 
produces 20 MeV of observable energy originates from a neutrino with larger energy, but just of 5 MeV on average. 
If the number of ES events in KII is larger, this could become more important. 



D. Speculations 

It is interesting to consider at this point some speculative scenarios, to investigate the question under which 
conditions we can increase the expected number of ES events: 




FIG. 5: Expected energy distribution of scattering events in KII, for £.j, = (0.5, 1, 2) x £e- Also shown is the probability to have 
one event with energy larger than 20 MeV. 
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• A distinguished astrophysical possibility is that there are main departures from a 'standard' collapse, and a 
large part of the emitted energy is not seen by inverse beta decay. Let us assume as an extreme case that the 
electron neutrinos have an average energy of 40 MeV and carry an energy of 1.5 x 10^^ erg The calculation 
reveals that the increase is not much larger: In KII, we expect Ndir — 0.59 rather than Ndir = 0.39. The reason 
is simply that oscillations transform the v,, into and Uti and the interaction cross section of these neutrinos 
is smaller. 

• Another possibility is to study which adjustment of the 'standard' scenario goes in the right direction. In 
particular, one can suggest that the are more energetic than what we assumed. This does not help for the 
number of events, but helps a bit to explain the fact that the directional events are among the most energetic 
ones. 

• The uncertainties in oscillations provide us with another degree of freedom. It seems to us very difficult to avoid 
the occurrence of MSW oscillations completely, but if ^13 is very small, we could get P^e — 0.3. The result 
Pee = sin^ 6'i2 = 0.3 assumes that the mantle of the star at densities of about 10 gr/cc was not essentially modified 
by the pre-coUapse events. The opposite case seems unlikely, but one could get Pee = 1 — sin^2f?i2/2 ^ 0.6. 
However, this does not help to increase Ndir with the 'standard' scenario, and it is of limited use to invoke 
non-standard scenarios with energetic z^e's, since in this case the reaction with oxygen are also called into play, 
see and Another case arises if 6'i3 is 'large' when the neutrino mass spectrum is inverted, rather than 
'normal' as considered previously in the text. In fact, in this hypothesis the IBD events are due to the flux 
and the flux of Ve (important for ES events) is 0.3'I>e + 0.7<1>2:. Thus, we are interested to consider the possibility 
of a large £e to increase the number of ES events, or the possibility that Se > ^xi since in this way we reduce 
the number of IBD events more than the ES events. 

• A more drastic attitude is to abandon completely the 'standard' idea of the collapse. A specific suggestion made 
in |37j| is that a large amount of neutrino radiation comes from 7r+ — > l^^i^/i decay. The largest contribution 
to scattering events comes from the electron neutrinos, that, due to oscillations originate exactly from (in 
fact, due to the loop-induced difference of potential between r and /x neutrinos —> Vi and a Ve happens to 
be produced with probability Pf^^e = 0.7). The are monochromatic with energy 29.8 MeV. If the energy 
injected in f^, is 5 x 10^^ erg we get about 1 ES event in 1MB and 3 ES events in KII (with a few additional 
Oxygen events). Apart from the obvious objection that we need to produce 10^^ (!) pions, we are left with the 
problem to explain the main part of the signal (that in the standard interpretation is attributed to IBD). 

In summary, we see that there are several interesting possibilities and the fact that we do not have a definitive theory 
of the collapse motivates their consideration, even though our cursory investigation seems to suggest that it is not so 
easy to produce radical modifications of the 'standard' paradigm. 

IV. SUMMARY AND DISCUSSION 

We reanalyzed the neutrino signal of SN1987A in 1MB and KII detectors in the light of new facts. In particular, 
occurrence of three neutrino oscillations (as defined in Sect lIB|l implies that the observed Pe have a 30 % contribution 
from the original i/^ or 9^-, while the observed electron neutrinos Ve are practically purely or i>r at the production. 
Thus, the 'standard' picture of the neutrino fluxes implies that the inverse beta decay signal is mostly sensitive to 
originally produced z/g, whereas an important contribution to the directional signal comes from i'^.^t- 

The hypothesis that most of the events were due to inverse beta decay is in agreement with the data, although 
the presence of one or more directional events is suggested by the shape of the angular distribution of the events. 
Even combining the information on the angular distribution with the a priori expectation for the number of events, 
an interesting hint that Kamiokande II dataset includes some elastic scattering events does remain (especially if £x is 
relatively large). 

It is conceivable that one can improve the agreement between the angular distribution and the expected (small) 
number of elastic scattering events by considering non-standard scenarios for the collapse. In the cases we considered 
the obtained improvements are interesting, but not dramatic. 

Let us conclude coming back to the 'standard' picture, and recalling which are the likely values of the parameters 
of the collapse we obtained: 

• In a 'standard' picture, we estimated from the data that the average energy of De is about Eq = {Ejy^) ~ 14 MeV. 
This is corroborated in particular by the average energy of 1MB events and by the fact that more events are 
observed at KII. Other pieces of data give contrasting hints: the angular distributions would like to have an 
average energy as high as possible, the average energy at KII suggest instead a lower average energy. In the 
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'standard' picture, we interpret these features as due to fluctuations, possibly with the contribution of one or 
more directional events in the data sets. As for the theoretical impact of this result, we note that a low average 
energy suggests an effective thermalization of the emitted antineutrinos. 

• The energy emitted in the collapse is about £b ~ 2.5 x 10^^ erg, for a distance of 52 kpc. Interestingly, this 
value is not far from simple minded theoretical expectations. Assuming further long wavelength oscillations in 
mirror neutrinos as in (4^ half of the neutrinos become invisible and £b should be doubled. Unfortunately, the 
calculations of £b do not seem to be precise enough to tell (2 — 3) from (4 — 6) x 10^'^ erg at present. 

• From the hint for elastic scattering event (s) we have some preference toward a comparably larger value of 
£x > Sg- This is compatible with current expectations, but it is unclear whether a large amount of radiation 
(that does not produce 'neutrino heating' for the delayed scenario) can be easily reconciled with the occurrence 
of the explosion, especially if this happens during the accretion phase. 

• Finally, it should be noted that there are hints (see and ^3) from astronomy that the distance of the Large 
Magellanic Cloud traditionally used is overestimated. If the new value of I? = 40 kpc is adopted, the energy 
emitted in neutrinos that we estimated has to be reduced by a factor of (40/52)^, namely f ~ 1.5 x 10^^ erg. 
Having little energy at our disposal is unlikely to help the occurrence of the supernova. This leads us to believe 
that the old determination of the distance is the correct one (as a matter of fact, more recent works [S^ 
argue from astronomical considerations that this is the case). 

To summarize our findings, we conclude that the 'standard' picture of neutrino emission and oscillations is not 
contradicted by SN1987A, and even more, the observed properties of the collapse seem to meet expectations. We 
believe that there is wide space for deviations from this picture, not only in consideration of the limited statistics but 
also due to certain features of the observed signals. From the discussion (and also in view of other considerations 
[Slf). it is evident that there is a great interest in obtaining larger samples of elastic scattering events and also of 
events due to from the next galactic supernova. 
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